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We investigate the effects of an external magnetic field in the gap structure of a color supercon¬ 
ductor with three massless quark flavors. Using an effective theory with four-fermion interactions, 
inspired by one-gluon exchange, we show that the long-range component B of the external magnetic 
field that penetrates the color-flavor locked (CFL) phase modifies its gap structure, producing a new 
phase of lower symmetry. A main outcome of our study is that the B field tends to strengthen the 
gaps formed by Q-charged and Q-neutral quarks that coupled among themselves through tree-level 
vertices. These gaps are enhanced by the field-dependent density of states of the Q-charged quarks 
on the Fermi surface. Our considerations are relevant for the study of highly magnetized compact 
stars. 
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After the suggestion that three-flavor quark matter 
may actually be the ground state of strong interactions 
[lj , quark stars were postulated as possible astrophysical 
objects. It is also very likely that quark matter occupies 
the inner regions of neutron stars. Our present knowl¬ 
edge of QCD at high baryonic density indicates that this 
new state of matter might be in a color superconduct¬ 
ing phase (for reviews, see [2|). On the other hand, it 
is well-known [j] that strong magnetic fields, as large as 
B ~ 10 12 — 10 14 G, exist in the surface of neutron stars, 
while in magnetars they are in the range B ~ 10 14 — 10 15 
G, and perhaps as high as 10 16 G The physical upper 
limit to the total neutron star magnetic field, as arising 
from comparing the magnetic and gravitational energies, 
is of order B ~ 10 18 G .3]- If quark stars are self-bound 
rather than gravitational-bound objects this upper limit 
could go higher. In this Letter we investigate the effect of 
a strong magnetic field in color superconductivity, with 
the aim of further studying its possible astrophysical im¬ 
plications. 

We will start by considering three massless quarks. In 
this case, it is well established that the ground state of 
high-dense QCD corresponds to the CFL (color-flavor 
locked) phase 0. In this phase, quarks form spin-zero 
Cooper pairs in the color- antitriplet, flavor-antitriplet 
representation, thereby breaking the original SU (3) co ior x 
SU (3) l x SU {3)rxU(1)b symmetry to the diagonal sub¬ 
group SU (3) co ior-f- l+r- One can now ask how this sce¬ 
nario will change when a magnetic field is switched on. 
Would the external field affect the pairing phenomena? 
In a conventional electromagnetic superconductor, since 
Cooper pairs are electrically charged, the electromagnetic 
gauge invariance is spontaneously broken and the pho¬ 
ton acquires a mass that can screen a weak magnetic 
field: this is the Meissner effect. In spin-zero color su¬ 
perconductivity, although the color condensate has non¬ 
zero electric charge, there is a linear combination of the 


photon and a gluon that remains massless jUj. The un¬ 
broken 17(1) group is generated, in flavor-color space, by 
Q = Q x 1 — 1 x Q, where Q is the electromagnetic 
charge generator :6]. Thus a spin-zero color superconduc¬ 
tor may be penetrated by a long-range remnant “rotated 
magnetic” field B. In the 9-dinrensional flavor-color rep¬ 
resentation that we will use in this paper, the Q charges 
of the different quarks are 
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in units of the Q-charge of the electron e = e cos@, where 
0 is the mixing angle |7j (we set h = c = 1 henceforth). 

Although the interaction of an external magnetic field 
with dense quark matter has been investigated by several 
authors , they disregarded the effect of the penetrat¬ 
ing B field on the gap structure. However, the B field 
can change the gap structure and lead to a new supercon¬ 
ducting phase. To understand how this can occur notice 
that due to the coupling of the charged quarks with the 
external B field, the color-flavor space is augmented by 
the Q charge operator, and consequently the order pa¬ 
rameter of the CFL splits in new independent pieces. 

Based on the above considerations, and imposing that 
the condensates should retain the highest degree of sym¬ 
metry, we propose the following ansatz for the color- 
flavor structure of the order paramete£of three massless 
quarks in the presence of a magnetic B field: 

A + = kt(U 0 -N)Clo + ktun 0 + k+Nn 0 +ktU[n + +Q.-\ 

.( 2 ) 

with color-flavor matrices defined as: Uq = 6 l a 6 3 b , 
U = Si5 l b , N = (616} + 6%6?)6j}6? + (a^b,i^ j)-, with 
a,b and i,j denoting color and flavor indexes respec¬ 
tively. The matrices = diag(l, 1 , 0 , 1 , 1 , 0 , 0 , 0 , 1 ), 
Q+ = diag(0, 0, 0,0, 0,0,1,1, 0), and = 
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diag(0, 0 , 1 , 0 , 0 , 1 , 0 , 0 , 0 ) are Q-charge projectors 
with algebra for 77 , 77 ' = 0,+,—, and 

t n_j_ t n_ = 1. 

An applied magnetic field reduces the flavor symme¬ 
tries of QCD, as only the d and s quarks have equal 
electromagnetic charge. Thus, the order parameter 
© implies the following symmetry breaking pattern: 
SU( 3) color x SU{2) l x SU{2) r x U(1) b x £/<">( 1)a x 
U( l) e . m . - SU(2) co1ot+l+r x U(l) e m . The U^\l) A 
symmetry is connected with the current which is an 
anomaly free linear combination of s, d and u axial cur¬ 
rents [ 13 . The locked SU( 2) corresponds to the maximal 
unbroken symmetry, and as such it maximizes the con¬ 
densation energy. Notice that it commutes with f/(l) e m . 

Therefore, there are 13 broken generators, 8 of which 
become the longitudinal components of massive gauge 
bosons, and 5 remain as Goldstone bosons. One is as¬ 
sociated to the spontaneous breaking of baryon symme¬ 
try, one with the breaking of the anomaly free 
and the remaining 3 to the breaking of the chiral SU (2) 
group. This symmetry breaking pattern suggests that 
the new phase has quantitative and qualitative differ¬ 
ences with respect to the CFL phase. We will call it 
magnetic CFL (MCFL) phase. In particular, the MCFL 
phase will possess a distinctive low energy physics. 

To trace back the physical origin of the new struc¬ 
tures in 0 we should take into account that despite the 
Q-neutrality of all the condensates, they can be com¬ 
posed either by neutral or by charged quarks. Conden¬ 
sates formed by Q-charged quarks feel the field directly 
through the minimal coupling of the background field B 
with the quarks in the pair. A subset of the condensates 
formed by Q-neutral quarks, can feel the presence of the 
field through the tree-level vertices that couple them to 
charged quarks. The gaps A^/g = (k n =F k c )/2 are anti¬ 
symmetric /symmetric combinations of condensates com¬ 
posed by charged quarks and condensates formed by this 
kind of neutral quarks. The gaps Aa/s = (fci T^2)/2, on 
the other hand, are antisymmetric/symmetric combina¬ 
tions of condensates formed by neutral quarks that do not 
belong to the above subset. The only way the field can 
affect Aa/s is through the system of highly non-linear 
coupled gap equations. The CFL gap matrix is obtained 
when = A a/s- In principle, the symmetries of the 

problem allow for two extra independent symmetric gaps. 
But these are only due to subleading color symmetric in¬ 
teractions, and are formed by neutral quarks that are not 
coupled to charged quarks, so they belong to the same 
class as As. Thus, in a first approach to the problem, we 
will consider that those can as well be described by As. 

To study the MCFL phase we use a Nambu-Jona- 
Lasinio (NJL) four-fermion interaction abstracted from 
one-gluon exchange [3]. This simplified treatment, al¬ 
though disregards the effect of the .B-field on the gluon 
dynamics and assumes the same NJL couplings for the 


system with and without magnetic field, keeps the main 
attributes of the theory, providing the correct qualitative 
physics. We will postpone the study within QCD for the 
future. 

The NJL model is defined by two parameters, a cou¬ 
pling constant g and an ultraviolet cutoff A. The cutoff 
should be higher than the typical energy scales in the sys¬ 
tem, that is, the chemical potential /i and the magnetic 
energy \fdB. 

The gap equation of the NJL model in coordinate space 
reads 

A+ = i T\ 2 X A r y >1 s 2 i(x,y)^tj, X A S (4) (x - y) , (3) 


where Aa and 7 ^ are the Gell-Mann and Dirac matrices, 
respectively. For simplicity, we have omitted explicit fla¬ 
vor, color and spinor indexes in the equation. S 21 is the 
21-component of the quark propagators in the Nambu- 
Gorkov representation. 

The computation of the field-dependent quark propa¬ 
gators is laborious (details will be given elsewhere llll| J. 
but it can be managed with the use of the Ritus’ method, 
originally developed for charged fermions 13 and re¬ 
cently extended to charged vector fields |l3j|. In Ri¬ 
tus’ approach the diagonalization in momentum space 
of charged fermion Green’s functions in the presence of a 
background magnetic field is carried out using the eigen¬ 
function matrices E p {x). These are the wave functions 
of the asymptotic states of charged fermions in a uni¬ 
form magnetic field and play the role in the magnetized 
medium of the usual plane-wave (Fourier) functions e lpx 
at zero field. With the help of the E p (x) functions, we 
first compute the propagators in momentum space, and 
then transform to coordinate space adequately. Leav¬ 
ing aside the color-flavor structure, the neutral quark 
propagators are of the same type as for the CFL phase. 
The charged (positive/negative) quark propagators in the 
background of a B field that lies in the 5-axis are 

g(±) __ A (±)J 5 kc __ A( ± )75 k c _ 

21 “ Pl-{\v (±) \-n?-kl Pl-{\v (±) \ + g) 2 ~kl 

(4) 


where = (0, ±\J2\eB\l,p3) are the spatial compo¬ 
nents of the momentum for (positive/negative) charged 
quarks and the integer number l labels the Landau levels. 

In 0, A±"' = (1 ± 707 • p l± ^)/2 are the energy projec¬ 
tors in the ultrarelativistic limit for (positive/negative) 


charged quarks in the external field, with p^ ' represent¬ 
ing the normalized charged-quark three-momentum. On 
the Fermi surface, the highest occupied Landau level is 

2 “I 

obtained as l max = „fi~, , where the bracket denotes 


21 eB\ 


integer part. 

Since color superconductivity, although of type I, al¬ 
lows the penetration of a rotated magnetic fieklAt is nat¬ 
ural to expect that the condensates made up of Q-cliarged 
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quarks will be strengthened by the nonzero B 1 because 
these paired quarks have opposite Q charges and opposite 
spins, hence parallel (instead of antiparallel) magnetic 
moments. The situation here has some resemblance to 
the magnetic catalysis (MC) of chiral symmetry break¬ 
ing 0 , in the sense that the magnetic field strengthens 
the pair formation. Despite this similarity, the way the 
field influences the pairing mechanism in the two cases is 
quite different. The particles participating in the chiral 
condensate are near the surface of the Dirac sea. The 
effect of a magnetic field there is to effectively reduce 
the dimension of the particles at the lowest Landau level 
(LLL), which in turn strengthens their effective coupling, 
catalyzing the chiral condensate. Color superconductiv¬ 
ity, on the other hand, involves quarks near the Fermi 
surface, with a pairing dynamics that is already (1 + 1)- 
dimensional. Therefore, the B field does not yield further 
dimensional reduction of the pairing dynamics near the 
Fermi surface and hence the LLL does not have a spe¬ 
cial significance here. Nevertheless, the field increases 
the density of states of the Q-charged quarks, and it is 
through this effect that the pairing of the charged parti¬ 
cles is reinforced by the penetrating magnetic field. Be¬ 
low we will analytically show that this is indeed the case. 

To solve the gap equation m for the whole range of 
magnetic-field strengths we need to use numerical meth¬ 
ods. We have found, however, a situation where an ana¬ 
lytical solution is possible. This corresponds to the case 
eB > y? /2. Taking into account that the leading con¬ 
tribution to the gap solution comes from quark energies 
near the Fermi level, it follows that for fields in this range 
only the LLL (/ = 0) contributes. 

Using the approximation Ag, Aa, Aa Ag, 

the gap equations decouple and the equation for A a is 

9 2 [ d 3 g A| 

3 A 2 A ( 2?r ) 3 - y) 2 + 2 (A ^) 2 

_i_ r j± _ a a 

3A2 j - ( 2 -) 2 ^ q -,r + (A*y ’ (j 

where the first/second term in the r.h.s. of Eq. (0 cor¬ 
responds to the contribution of Q- neutral/charged quark 
propagators, respectively. For the last one, we dropped 
all Landau levels but the lowest, as we are interested in 
the leading term. 

The solution of Eq. 0 reads 

3A 2 t r 2 

g 2 (^y 2 + eB 

with 6 = A — y, to be compared with the antisymmetric 
CFL gap @ 

A^~2,/^exp(-|tiL). ( 7 ) 


A^ ~ 2 \fdy exp ^ 



In this approximation the remaining gap equations 
read 
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( 8 ) 


(9) 


( 10 ) 


We express below the solution of these gap equations as 
ratios over the CFL antisymmetric and symmetric gaps 


Aa 1 / 36 21 1 

A^ FL 2U/ 34 ) e P V 17^; 17 2 j(l -|- y'j 

where x = g 2 y 2 / K 2 tt 2 ^ and y = eB/y 2 , and 

Af Ag/3 9 y -l \ 

a c fl A Cfl W t 2xln2y + l) 

A s Aa 3/ 4 \ 

a c F l A cfl 2 V i + yj • 



( 11 ) 


( 12 ) 

(13) 


Note that our analytic solutions are only valid at strong 
magnetic fields. The lower value eB ~ y 2 /2 corresponds 
to eB ~ (0.8 - 1.1) • 10 18 G, for y ~ 350 - 400 MeV. 
For fields of this order and larger the A^ gap is larger 
than A^fl at the same density values. Note also that for 
eB > A* 2 , A^, g grow and A A / S tend to lower, and they 
clearly split. How fast or slowly they do depends very 
much on the values of the NJL couplings. For example, 
for x ~ 0.3 0, one finds Aa ~ 0.2 A a for y = 3/2, 
while for x ~ 1 then Aa ~ 0.5 A a . 

All the gaps feel the presence of the external magnetic 
field. As expected, the effect of the magnetic field in A® 
is to increase the density of states, which enters in the 
argument of the exponential as typical of a BCS solution. 
The density of states appearing in © is just the sum of 
those of neutral and charged particles participating in the 
given gap equation (for each Landau level, the density of 
states around the Fermi surface for a charged quark is 
eB/2n 2 ). 
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All the Q-charged quarks have common gap A((. 
Hence, the densities of the charged quarks are all equal. 
As two of these quarks have positive Q charge, while 
the other two have it negative, the Q neutrality of the 
medium is guaranteed without having to introduce any 
electron density. 

Our zero-temperature results imply that a propagat¬ 
ing rotated photon with energy less than the light¬ 
est charged quark mode cannot scatter, since all the 
Q-charged quarks acquire a gap and all the Nambu- 
Goldstone bosons are neutral. The anisotropy present 
in the background of an external magnetic field and the 
existence of charged Goldstone bosons in CFL but not 
in MCFL indicates a rather different low energy physics, 
including transport properties. In particular, the MCFL 
superconductor is transparent and behaves at T = 0 as an 
anisotropic dielectric, as opposed to the isotropic dielec¬ 
tric behavior of the CFL phase mu However, similar 
to the CFL, the medium will become optically opaque as 
soon as leptons are thermally excited [lZ(. 

In previous analysis (see first paper in Sj) the criti¬ 
cal magnetic field at which the CFL pairing is destroyed 
was estimated to be ~ 10 20 G. This estimate was based 
in a field-independent CFL pairing energy 
Considering that in the MCFL phase A® increases with 
the field, it is natural to expect that the critical field will 
be even larger than 10 2O G. Since such extremely strong 
fields will surpass all the energy scales of the system, the 
quark infrared dynamics will become predominant, and 
the phenomenon of magnetic catalysis of chiral symme¬ 
try breaking 13 will be activated, producing a phase 
with quark-antiquark condensates but no quark-quark 
condensate. These two phases will have to be connected 
by a phase transition, as they have different number of 
Goldstone bosons due to the breaking of baryon symme¬ 
try, which only occurs in the superconducting phase. 

Let us stress that in this work we have not considered 
the implications of finite quark masses. This, together 
with a careful study of the effects of the magnetic field 
in the low energy physics, in transport properties or in 
neutrino dynamics will be the subject of future investi¬ 
gations. 

In conclusion, we have found that a magnetic field leads 
to the formation of a new color-flavor locking phase, char¬ 
acterized by a smaller vector symmetry than the CFL 
phase. The essential role of the penetrating magnetic 
field is to modify the density of states of charged quarks 
on the Fermi surface. To better understand the relevance 
of this new phase in astrophysics we need to explore the 
region of moderately strong magnetic fields eB < /i 2 /2, 
which requires to carry out a numerical study of the gap 
equations including the effect of higher Landau levels. 
Because the total density of states around the Fermi sur¬ 
face for charged particles does not vary monotonically 
with the number of Landau levels, we still expect to find a 
meaningful splitting of the gaps at these fields and there¬ 


fore a qualitative separation between the CFL and MCFL 
phases. 
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